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t-H . 

f^**) ■ We present measurements of absorption by the 21cm hyperfine transition of neutral hydrogen 

toward radio sources at substantial look-back times. These data are used in combination with 
observations of rotational transitions of common interstellar molecules, to set limits on the evolution 
of the fine structure constant: ^ < 3.5 x 10 -15 year -1 , to a look-back time of 4.8 Gyr. The 
neutral hydrogen observations employed Very Long Baseline Interferometry in order to mitigate 
the substantial uncertainty arising from the fact that observations at very different wavelengths 

qq ' may probe different lines-of-site due to frequency dependent structure of the background source. 

We discuss the implication of these results on theories unifying natural forces based on compact 

£f~) . quantum dimensions. In the context of string theory, the limit on the secular evolution of the 

scale factor of these compact dimensions, R, is f < 10~ 15 year . Including terrestrial and other 
astronomical measurements places limits (2a) on slow oscillations of R from the present to the epoch 
of cosmic nucleo-synthesis, just seconds after the big bang, of < 10 5 . 
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Since the startling realization that we live in an evolving universe, there have been numerous hypotheses concerning 
the cosmic evolution of physical constants. Dirac pointed out that the ratio between the strength of the gravitational 
force to that of the electromagnetic force « 10 39 ~ the age of the universe measured in atomic units of time. (ie. 'le 
tempon' = light crossing time of the classical radius of the electron). He speculated that such large dimensionless 
numbers may be more fundamental than the supposed fundamental constants, and that the constants evolve in 
such a way that observers at any given epoch reach similar conclusions. This line of reasoning leads to a variation 
of the gravitational constant]], G, over cosmic time of order: ^ ~ H Q ~ 10~ 10 year -1 , where H D is the Hubble 
constant quantifying the local expansion of the universe. A similar prediction on G comes from the Brans-Dicke 
equation relating mass/energy to a cosmic scalar field ||. Variation of G at the Hubble-rate has been ruled out by 
measurements of planetary orbits Gamow || speculated that the charge on the electron, e, or more precisely, 

the fine structure constant, a, varies and not G. This would lead to variation of the Rydberg constant and hence 
contribute to the observed redshift, z, of spectral features in distant objects normally attributed to the evolution of 
the 4-dimensional (4D) cosmic scale factor. Lastly, Sisterna and Vucetich || point out that stringent limits to the 
variation of physical constants support Einstein's strong equivalence principle, and hence provide 'accurate verification' 
for general relativity as the correct low-energy theory of gravity. 

The idea of cosmic variation of physical constants has been revisited with the advent of models unifying the forces 
of nature based on the symmetry properties of quantum dimensions, such as the Kaluza-Klein (KK) model, or the 
more general requirement of extra dimensions in superstring theory (SS) These extra dimensions have a scale 

factor, R, of order the Planck scale, R ss 10 -33 cm, and manifest themselves only during the first instant of creation, 
corresponding to the Planck time, 10 -43 seconds after the big bang, or at energies above 10 19 GeV. These compact 
dimensions quickly vanish during the cosmic expansion of our familiar 4D space-time, however they may still have 
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1 Study of the variation of fundamental constants must consider dimensionless forms involving products of the constant in 

Gm 2 

question with h, m p , and c, since dimensionless quantities are invariant under coordinate transformation, for example: 7 = Rc p 
or a = j— Bpl. We use H = 75 km s" 1 Mpc" 1 , and a deceleration parameter, q = 0.5. 
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observable consequences, since the constants of nature observed in 4D are the result of integration over the extra 
dimensions It has been hypothesized that a variation of R with cosmic epoch could lead to a variation of the 
physical constants measured in 4D PJlQfl. In particular, the fine structure constant, a, is predicted to behave as 
RT 2 in KK theories, and R~ 6 in SS theories p||ll[|. Unfortunately, while the time dependence of the cosmic scale 
factor in 4D space-time is well understood in terms of Einstein's theory of gravity, there is no analogous prediction 
concerning the evolution of R in extra-dimensional theories. It is possible that R increases or decreases monotonically, 
or even oscillates, with cosmic time JTo| . Still, a measurement of the cosmic evolution of a would provide qualitative 
supporting evidence for the existence of compact dimensions fl2|| . 

Limits to the evolution of a include laboratory measurements ]13| ] , consideration of the abundances of radioactive 
isotopes jl3| , and consideration of fluctuations in the microwave background and other cosmological constraints Jl4|,[l5| ■ 
The most stringent terrestrial measurement comes from the Oklo natural reactor, which occurred about 1.8 Gyr ago. 
The 95% confidence limit (2a) from these calculations is |~| < 7 x 10~ 17 year -1 p7| . Calculations of primordial 
nucleo-synthesis allow for a maximum variation of |^| < 11 x 10~ 15 year _1 ( JTT|,p|| although see [^6|). Accurate 
spectroscopy of absorption and emission lines from objects at cosmologically significant redshifts can be used to set 
limits on the evolution of physical constants Jl8| p0| , ^2|j5[| , Such studies lead to limits of p2] , pl[ : |^| < 2 x 10~ 15 year -1 . 

Drinkwater et al. 1 23 1 and Wiklind and Combes p4| have compared absorption by molecular rotational transitions 
at millimeter wavelengths to HI 21cm absorption to determine the evolution of the product Y = g p a 2 , where g p is 
the proton-to-electron magnetic moment. These cm and mm measurements have the potential advantage over optical 
spectroscopy in that spectral resolutions of 1 km s _1 or better are easily obtained, and the absorption lines themselves 
can be extremely narrow |24f| , with the ultimate limitation 10 -8 ) being the accuracy of the lab measurements of the 
transitions in question. Based on errors due to signal-to-noise alone, the current best limit (2<r) using this technique 
is |^| < 1 x 1CP 15 year -1 . However, these measurements have a potentially much larger uncertainty arising from 
possible differences in the velocities of the HI 21cm and molecular absorbing gas within a given galaxy ]2^ , |l2| , pil"| . 
Velocity differences can arise both along a given line-of-site, and also due to the fact that observations at very different 
wavelengths (eg. mm versus cm) may probe different lines-of-site due to frequency dependent spatial structure of the 
background source. If line-of-sight differences occur on kilo-parsec (kpc) scales, then systematic velocity differences 
can arise due to the galaxian potential, and can be of order 100 km s _1 . If line-on-site differences can be limited to 
sub-kpc scales, then it can be argued that the residual uncertainty is likely to be of order 10 km s -1 , ie. comparable 
to the typical velocity dispersion of the interstellar medium in galaxies |25j] . 

In this paper we present observations of HI 21cm absorption at intermediate redshifts using the technique of Very 
Long Baseline Interferometry (VLBI) |^7j . These observations provide spatial resolutions of tens of milli-arcseconds 
(mas), corresponding to sub-kpc spatial scales, and hence mitigate the potential problem of probing different lines- 
of-sight at different wavelengths. In the following analysis we will quote heliocentric redshifts, z, and differences in 
redshifts, Az. In terms of observed frequency differences, Ais , source-frame velocity differences, Av, and the variation 
of the physical quantity in question, AY: 4^ = — ~ 

We have observed the HI 21cm absorption toward the cosmic radio sources 0218+357 and 1413+135 using VLBI 
techniques (Fig. 1). The source 0218+357 is a gravitationally lensed background radio source, and the absorption is by 
gas in the lensing galaxy ]2S|j2"g| ] . The source 1413+135 is a radio loud AGN at the center of an edge-on spiral galaxy, 
and the absorption is by gas in a molecular cloud in the disk of the parent galaxy 30 24j] . The data were processed 
at the VLBA correlator in Socorro NM, and images were generated using the Astronomical Image Processing System 
(AIPS). The data were fringe fit starting with a point source model for 1413+135 and a double source for 0218+357, 
and then imaged using standard hybrid imaging techniques |27j| . In each case spectral data image cubes of 256 
channels over 4 MHz bandwidth were analyzed. The details of the observations and their astrophysical implications 
will be presented elsewhere. Herein we concentrate on the information relevant to constraining the evolution of a. 

For 0218+357 the continuum image has a resolution of 80 mas, corresponding to a physical scale of 390 pc, and shows 
two lensed radio components separated by 350mas. We find that the dominant HI 21cm absorption line is toward the 
southwest component. This is also true for the molecular absorption plfl . We have fit a two component Gaussian 
model to the absorption spectrum using a least squares algorithm. The principal component has a peak optical depth, 
t = 0.12, and a Full Width at Half Maximum, FWHM = 15 km s _1 , at a redshift of z = 0.684676 ± 0.000005. The 
error corresponds to a change in z leading to a unit change in reduced x 2 ■ The peak of the molecular absorption in 
this system is at z — 0.684693 ± 0.000001 pij] , hence the difference between the molecular and HI 21cm redshifts is: 
4^- = 1.0 x 10 -5 ±3 x 10 -6 . 

1 + 2 

For 1413+135 the continuum image has a resolution of 20 mas, corresponding to a spatial scale of 70 pc, and shows 
a core-jet morphology extending over 60 mas. The position of the inverted-spectrum nucleus is shown as a large cross, 
and again, at mm wavelengths only the nuclear radio component is detected [ p2[ . At 1143 MHz the jet dominates 
the total radio continuum emission from 1413+135. The absorption spectrum at the peak surface brightness of the 
jet, indicated by a small cross, is shown in Figure 1. The sensitivity of these observations is insufficient to detect 
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absorption toward the weaker nuclear radio component. We fit the spectrum of the jet with a one component Gaussian 
model, resulting in r = 0.89 and FWHM = 15 km s" 1 , at z = 0.246693 ± 0.000001. The molecular absorption in this 
system is extremely narrow (FWHM w 1 km s -1 ) with a peak at z = 0.2467091 ± 0.0000003 Q, hence the difference 
between the molecular and HI 21cm redshifts is: = 1.3 x 10~ 5 ± 8 x 10~ 7 . 

The errors given above are la strictly based on the fitting process and the signal-to-noise ratio (SNR) in the observed 
spectra. There are a number of other systematic errors in the measurement and Gaussian fitting process which must 
be considered, including: (i) the transfer of the sky frequencies from the telescopes through the correlator to the data 
reduction programs, and (ii) the interpolation process when deriving redshifts to accuracies better than the spectral 
channel width of 15 kHz, including possible variations in the zero intensity level, the number of Gaussians used in the 
modeling, and the assumption of a Gaussian shape for the line profile. The results indicate that the measurement 
errors other than SNR are < 10 kHz, corresponding to 3 km s — 1 , or < 1 x 10~ 5 . 

The sum total of the uncertainties discussed above are comparable to the most accurate astronomical measurement 
to date p^ , pl| , and could be improved with higher SNR, higher spectral resolution observations. However, there 
remains the possibility of systematic velocity offsets between the HI 21cm absorbing gas and the molecular absorbing 
gas in a given galaxy. Drinkwater et al. |23| demonstrate that the velocity of HI 21cm absorption in the Galaxy 
shows a statistical correlation with that of molecular absorption, with a dispersion of only 1.2 km s -1 . On the other 
hand, there are clearly Galactic and extragalactic examples in which significant velocity differences occur [^4|, and 
a statistical analysis of the type in Ref. |23| will be rigorously valid only when applied to a statistical sample of 
absorbers. Note that this error does not apply to calculations of — based on fine structure absorption doublets by 
alkaline metals [p2|. 

For 1413+135 the molecular and HI 21cm absorption spectra are narrow in velocity, suggesting absorption by the 
general ISM of the host galaxy. But the HI 21cm absorption occurs along a line-of-site toward the jet which is offset 
from that toward the nucleus by 25mas. This line-of-site offset will give rise to a velocity difference less than 10 km 
s -1 as long as the absorbing gas is beyond 0.5 kpc from the nucleus, assuming a flat rotation curve for the galaxy 
with a velocity of 200 km s _1 . For 0218+357 the absorption is again narrow, and our VLBI observations probe 
a region of diameter < 390 pc through a face-on spiral galaxy p8| . In summary, it is reasonable to assume that 
the dominant uncertainty between the molecular and HI 21cm redshifts in both 1413+135 and 0218+357 is due to 
small-scale (sub-kpc) ISM motions, i.e. ss 10 km s -1 . A velocity uncertainty of 10 km s _1 leads to a limit to the 
evolution of the fine structure constant of | — | < 3.5 x 10~ 15 year -1 to a look-back time, t, of 4.8 Gyr for 0218+357, 
and | — | < 6.7 x 10~ 15 year -1 to t — 2.5 Gyr for 1413+135, assuming that g p is constant |Q. These limits are 
sufficient to rule-out hypotheses D and E summarized in Dyson pj, in which a oc t, and a cx log(t), respectively. 

In Table 1 we summarize the limits to — over time, and to -§ in KK and SS theories. The most stringent limit to 
secular evolution of R remains the Oklo limit, although this limit has been called into question recently |9j. Hence, 
limits using other methods provide an important check of the Oklo limit. Also, all the limits taken together argue 
strongly against a slowly oscillating R with < 10 -5 over the entire history of the universe. And the limits 
to — all have different functional dependencies on a and other physical constants. The lack of variation observed 
for any of the different products of physical constants argues against models of a 'cosmic conspiracy' in which the 
individual constants vary in concert to result in a given observable remaining invariant. One method for overcoming 
the uncertainty due to small scale ISM motions is to observe a large sample of sources (> 100), and rely on the 



statistical correlation described in Drinkwater et al. |23 . Study of such a large sample will only become possible with 



the greatly increased sensitivity of the Atacama Large Millmeter Array, and the Square Kilometer Array. 
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TABLE I. Limits (2 a) to the evolution of the fine structure constant and the scale factor of compact dimensions. Quantities 
with redshifts are based on astronomical spectroscopy. 



Method 


z 


Look-back Time 
Gyr 


xicr 6 


1 ( ss ) 

xl0~ 16 year" 1 


| (KK) 
xlO" 16 year" 1 


Reference 


Laboratory 




10" y 


4 x 10" s 


70 


200 




13 




Oklo Reactor 


- 


1.8 


0.13 


0.12 


0.35 




17 




1413+135 


0.25 


2.5 


17 


11 


33 


This paper 


0218+357 


0.68 


4.8 


17 


6 


18 


This paper 


Radioactive elements 




5 


24 


8 


25 


1 


11 




Fine Structure 


0.8 


5.2 


8 


2.6 


7.7 


1 


22 




Fine Structure 


1.3 


6.3 


8 


2.1 


6.3 


1 






Metal + HI 21cm 


1.78 


7.0 


15 


3.5 


10 


1 


21 




Primordial Nucleosynthesis 




8.9 


100 


18 


55 




11 





FIG. 1. The contour images (left) show total radio continuum surface brightness distributions for 0218+357 (top) and 
1413+135 (bottom) derived from VLBI observations at the frequency of the redshifted HI 21cm absorption lines. The observa- 
tions of 0218+357 were made using the European VLBI Network in October 1998, while those for 1413+135 were made using 
the Very Long Baseline Array operated by the National Radio Astronomy Observatory during July 1998. The spatial resolu- 
tions are shown in the insets (FWHM). The total flux density of 1413+135 at this frequency is 1.2 Jy while that for 0218+357 
is 1.8 Jy. The contour levels are logarithmic in the square root two with an arbitrary absolute scale since no a priori gain 
calibration was applied due to lack of accurate system temperature measurements at these low frequencies. The spectra (right) 
show optical depth in redshifted HI 21cm absorption toward the positions designated by the crosses in the continuum images. 
The solid line is the measured data, and the long dash line is the Gaussian model fit. In the case of 1413+135, the position of 
the AGN is shown with the large cross, but the HI 21cm absorption spectrum was made at the position of the peak surface 
brightness of the jet, corresponding to the smaller cross. The frequency scale has been corrected to a heliocentric rest-frame. 
The short dash lines show the CO absorption line profiles, scaled in frequency accordingly and on an arbitrary optical depth 
scale, as observed by the IRAM 30m telescope (see [25] and references therein). The arrows indicate the expected HI frequencies 
based on Gaussian model fits to the molecular absorption lines. For 0218+357 the profile displayed is for CO(l-2), while that 
for 1413+135 is CO(0-1). 
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